MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
ARTIFICIAL INTELLIGENCE LABORATORY 


A.I. Memo No. 657 


September, 1981 


Nonlinear interactions in a dendritic tree: 
localization, timing and role in information processing 

T. Poggio and C. Kocli* 

*Max-Planck-Institut fur biologische Kybcrnctik 
Spcmannstrasse 38, D 7400 Tubingen, FRG 


Abstract: In a dendritic tree transient synaptic inputs activating ionic conductances with an equi¬ 
librium potential near the resting potential can veto very effectively other excitatory inputs. Analog 
operations of this type can be very specific with respect to relative locations of the inputs and their 
timing. We examine with computer experiments the precise conditions underlying this effect in the 
case of a (5-like cat retinal ganglion cell. The critical condition required for strong and specific interac¬ 


tions is that the peak inhibitory conductance change must be sufficiently large, i. e. about 5 • 10“ 8 5 
or more, almost independently of other electrical parameters. In diis case, a passive dendritic tree 
may perform hundreds of independent analog operations on its synaptic inputs, without requiring any 
threshold mechanism. 
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Figure 1 . A cal retinal ganglion cell of the 6 type (Boycott and Wassle, 1974, Poggio ct al., 1981; Koch el al, 
1982). In the calculations reported in Figs. 2 and 3, C m = 2 nFcm~ 2 , R m = 2500ficm 2 , IR = 700 cm. 


When two regions of dendritic membrane in close electrical proximity experience simultaneous 
conductance changes • induced by synaptic inputs - the resulting postsynaptic potential at the soma is 
not the linear sum of the potentials generated by each synapse alone. The existence of such nonlinear 
interactions in a passive dendritic tree has been long recognized, both theoretically and experimen¬ 
tally (sec Rinzel and Rail 1974; Redman 1976; Barrett and Crill, 1974). Recently it has been proposed 
that this may represent an important mechanism for implementing elementary information processing 
operations locally in a passive dendritic tree (Poggio and Torre, 1978,1981). In particular, Torre and 
Poggio (1978,1981) have suggested that direction selectivity of some retinal ganglion cells is due to the 
nonlinear interaction on the dendritic membrane of the cell of an excitatory conductance change and 
of an inhibitory input with an equilibrium potential near the resting potential (shunting inhibition). 
We have analyzed the interaction of a transient inhibitory conductance input of the shunting type 
with an excitatory input, on the basis of cable theory, for various types of cat retinal ganglion cells. 
In this note we report, in particular, our analysis for the 6 cell shown in Fig.l, whose quantitative 
morphology was obtained from histological material of Boycott and Wassle (1974). We wish to show 
that inputs of the shunting type can “veto” in a nonlinear way die depolarising effect on the soma of 
an excitatory input. 


This nonlinear interaction can be 
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(a) strong 


(b) specific with respect to the relative position of excitation and inhibition 


(c) tuned to their (relative) timing. 

Properties (a) to (c) depend on the electrical parameters and on the morphology of the cell. By far 
the most critical condition is the size cf the conductance changes involved. For our cells properties 
(a) to (c) arc certainly valid if the peak conductance change is in the order of 5 • lO -8 ^ or larger, 
for a wide range of values of the other electrical parameters. Among die various types of cat retinal 
ganglion cells (a, /?, 7 and S according to Boycott's and Wacsslc’s classification (1974)) this inhibition 
is maximal for 7 and S cells, while it is relatively weaker for the a and (3 cells (Koch et a 1 . 1982). 
The branching structure, die length and the diameters of each dendritic segment were determined 
from histological material of Boycott and Waessle, as described in detail (Poggio et al., 1981; Koch 
et al., 1982). The dendritic tree was approximated in terms of short segments, each being equivalent 
to a cylinder. An algorithm (modified from Butz and Cowan (1974)) was implemented in a program 
to compute from these data (for a range of values of the electrical parameters, ( C rn ,R m ,Ri )) the 
linear electrical properties of die cell. The main assumption is that die membrane is passive (and diat 
linear cable theory holds). The program computes the complex transfer resistances for any 

two locations i,j in die dendritic tree. If a current Ij is injected at location j die resulting voltage at 
location i is given by 


vm=m * wi) (i) 

where * indicates convolution and Kij(t) is the inverse Former transform of Kij(u>). The set 
of for various locations i,j characterizes completely die (linear) electrical properties of a 

branched passive cable. We made use of two general properties of the Kij(u >) for dendritic trees (cf. 
Koch et al. 1982): 

symmetry, i.e. 

kijW = kji(u>) ( 2 ) 


separability, i.e., 



kn{u>) 



if l is on the “direct path” from i to j. 

Let us consider die case of an excitatory synapse modulating the conductance g e to an ionic species 
with equilibrium potential E e > V rest in location e and an inhibitory synapse modulating the conduc¬ 
tance gi to an ionic species with equilibrium potential £) = V rest — 0 in location i (for evidence of 
shundng inhibition in ganglion cells see Baylor and Fettiplace 1979). For transient conductance inputs 
the system of Volterra integral equations giving the resulting somatic potential is 
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Figure 2. The depolarization in the soma of the delta cell of Fig. 1 for an excitatory input g e at location 

—t 

1 and an inhibitor}' input g t at location 8. Both inputs have the time course g{i) = t 4 e t P cak and are shown 
in fig. 2a. Location and timing of inhibition are optimal (i.e. the inhibition is delayed by 1.5mscc see Fig. 
3). The excitatory battery is E c = 80 mV and the inhibitory battery is = OmV (relative to the resting 
potential). For both inputs, time to peak is t p = 25msec, g max — 10” 7 S\ Fig. 2b shows the corresponding 
somatic depolarization in the absence (\'e) and in the presence (Ve+i) of inhibition. The maximum of V in 
Fig. 2b is for t = 31 msec. Since g(t) peaks at t = 25 msec, the traveling time from location 1 to the 
soma is about 6 msec, which equals the phase-time-lag of the transfer function K is (u>) for cu = 0. Inhibition 
alone is “invisible” (because E{ = 0); its effect appears only when simultaneous excitation takes place, as 
expected for a nonlinear interaction. 


v.{t) = g e (t)[E e - 

- V e (t)} * K es {t) - 

(9i(t)Vi(t)) * K is (t) 

V e (t) = g t {t)% - 

- V e (<)l * K ee (t) - 

G 9i(t)Vi(t )) * K ie {t) 

V i {t) = g e {t)[E e - 

- V e (t)] * K ei (t) - 

mm) * K i&) 


* 

(4) 


This system of equations was integrated numerically for given inputs g e and gi using die Kij func¬ 
tions calculated by our program for the specific cell. Figure 2 shows the somatic potential V 8 for 
various inputs g e and p t . 
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Delay zltCmsec] 


Figure 3. F factor (ratio of the maximum of the somatic depolarisation without inhibition to the somatic 
depolarisation with inhibition) for various locations of excitation and (shunting) inhibition in the cell of Fig. 
1 as a function of relative timing. Ihc conductance inputs are as in Fig. 2, but with various delays between 
them. Notice that the optimal timing for e = i is not for At = 0 (Segev,pcrs. comm.; Poggio and Torre, 
1981); the deviation is, however, very small. In these calculations, we have assumed C 7n = 2jiFcro— 2 , 
Ri — 700cm, R m = 25000cm 2 . The excitation is always at location 1, while inhibition can be behind 
excitation (7), coincident with excitation, on 'the direct path' (8)or at the soma ( S ). The basic results hold 
over a wide range of parameters; in particular R m can be increased by several order of magnitudes without a 
significant change in the “direct path" property, which depends critically on branching geometry and effective 
intracellular resistance. 


A simple measure of the effectiveness of shunting inhibition is the ratio (F) between the maximum 
of somatic depolarization in the absence of inhibition and in the presence of the inhibitory input. 
Equations (4) can be solved analytically for steady state inputs. In this case, it can be shown, using the 
properties above that 



9e^es 

1 9eFee 


1 -(- g<K ee -f- 9 jKu 4- g t gjK* 
g e K ea + g e giK+ 



with K* = K ee Ku — KieKie and K+ — K es K li — K ie K ia . 

(The Kij represent the DC values of the transfer impedances; that is 0)). 

This formula shows that F can take quite high values, even for & significantly smaller than g e , quite 
differently from the lumped circuit case (see Koch et al, 1982). Figs. 2 and 3 (for tire delta cell) 
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show that sufficicnly large inhibition is highly effective in vetoing excitation also for transient inputs, 
provided location and timing are optimal. For steady state inputs it is possible to prove rigorously 
(from cq. 5) that the most effective location for inhibition is always on tire direct path from the 
location of the excitatory synapse to the soma (optimal location on the path depends on the values of 
ft, and ft, approaching e for small g c , ft and s for large ft). For ft > 5 • 10~~ 8 S, F is quite large, for 
"on path" inhibition over a wide range of R m (between 500ficm 2 and 1 MQcni 2 ). F increases for 
increasing R rn but not very much (at most by a factor 3 for R m increasing 3 orders of magnitude). 
Over tire same R, n range the "on path" effect maintains a good specificity. The situation changes for 
conductance peak values below 5 • 10~ s S. 


Two cases must be then distinguished: (a) R m is low: then F is between 1 and 2 (foi R m = 
500ncm 2 , ft = 10“ 9 S' and ft, = 1(P 8 SF = 1.6 under optimal conditions), (b )R m is larger: then 
the F values arc also larger (for R m = 20000ficm 2 and tire above conductance values F == 2.9) 
but the "on the path" effect is weak. The effect of inhibition depends in this case mainly on the 
distance from excitation (especially for small values of ft ); for very large values of R rn (around IMfl 
) inhibitory inputs at any location throughout the dendritic tree have very similar effect. Notice that 
for "large" values of R rn the soma usually becomes the optimal location for inhibition. F values and 
the "on path" specificity arc more sensitive to changes in R,. increasing with intracellular resistance. 
The physiological range of R{ is however quite restricted (between 50 and lOOQcm). Thus one can 
envisage that a neuron can work in two different modes of operation, depending on the strenght of the 
synaptic inputs. For small conductance inputs their interaction is mainly a function of the distance be¬ 
tween the synapses involved, while for larger inputs the interaction is more specific, showing a strong 
"on path” effect. Fig. 2b shows an example of the strength and specificity of shunting inhibition for 
the transient inputs shown in Fig. 2a. For these conductance changes inhibition on the direct path is 
strong whereas inhibition behind excitation or on a side branch a few microns further away is almost 
completely ineffective. Fig 3 shows how tire timing of die excitatory vs. the inhibitory input influences 
the effectiveness of the interaction. The width of the tuning curves of Pig. 3 reflects the time course 
of tire input conductance changes (shown in Fig. 2) and the cable properties. The optimal delay 
is essentially due to tire propagation time from excitation to the location of inhibition (die effective 
propagation time for “long” transients is given by die phase-time-lag of F' i j(w)| u)==0 , i.e., by 



(see Koch and Poggio, in preparation). For t pca k — 25 msec, inhibition "on the path" can effectively 
veto excitation if it occurs within 10msec of die onset of die excitatory input. Similar effects can be 
obtained by a longer lasting inhibition (instead of delayed inhibition). We investigated die effect of 
changing t pea k of the inhibitory input for a fixed delay between excitation and inhibition. Reducing 
tpeak for inhibition below die fixed value for excitation reduces drastically F. Thus, inhibition 
in order to be effective must last at least as long as excitation, but does not need to last much 
longer (depending on location).When the locations of excitation and inhibition coincide, maximal 
effectiveness of inhibition is reached for a t pea k as long as the excitatory t pea k- Inhibition in the soma 
needs to last roughly twice as long as excitation, in order to be maximally effective. 

Because of the strength and specificity of such nonlinear interactions we wish to propose that 
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tlicy may perform characteristic information processing operations in (passive) dendritic trees. Since 
inhibition vetoes effectively more distal excitatory inputs only when it is “on the path” to the soma, 
a variety of local operations can be performed, exploiting the branching geometry of a dendritic 
tree with a suitable localization of excitatory and inhibitory inputs. Timing of inputs provides an 
additional important control variable: inhibition “on the path” can veto in an “and-not” fashion an 
excitatory input only when it takes place within a well defined temporal window. 


Our results hold also for the more unusual case of a conductance decreasing inhibitory input 
(Dudcl and Kuffler, 1960; Gcrschenfeld and Paupardin-Tritsch, 1974): in this case, an inhibitory 
input facilitates the excitatory effect implementing an analog approximation of an “and” logical 
operation instead of the “and-not” discussed in this paper. 

Simple operations of the "and not” type may underly, for instance, direction selectivity to motion 
of certain neurons (Torre and Poggio,1978; Poggio ct al., 1981; Koch et al., 1982); mechanisms of the 
"and" type could be used in other motion sensitive cells. We recently proposed (Poggio ct al., 1981; 
Koch et al., 1982) that a 5-likc morphology is tire substratum of directional selectivity in the cat retina. 


In summary, we have shown that shunting inhibition may veto ox facilitate very effectively and 
specifically excitatory inputs in a 5-like cell, if the associated conductance change is large enough. 
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